UG-FT-184/05, CAFPE-54/05, ROMA-1403-05 
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Cosmogenic neutrinos reach the Earth with energies around 10 9 GeV, and their interactions with 
matter will be measured in upcoming experiments (Auger, IceCube) . Models with extra dimensions 
and the fundamental scale at the TeV could imply signals in these experiments. In particular, the 
production of microscopic black holes by cosmogenic neutrinos has been extensively studied in the 
literature. Here we make a complete analysis of gravity-mediated interactions at larger distances, 
where they can be calculated in the eikonal approximation. In these processes a neutrino of energy 
E v interacts elastically with a parton inside a nucleon, loses a small fraction y of its energy, and starts 
a hadronic shower of energy yE v <C E v . We analyze the ultraviolet dependence and the relevance 
of graviton emission in these processes, and show that they are negligible. We also study the energy 
distribution of cosmogenic events in AMANDA and IceCube and the possibility of multiple-bang 
events. For any neutrino flux, the observation of an enhanced rate of neutral current events above 
100 TeV in neutrino telescopes could be explained by TeV-gravity interactions. The values of the 
fundamental scale of gravity that IceCube could reach are comparable to those to be explored at 
the LHC. 



PACS numbers: 04.50.+h, 13.15.+g, 96.40.Tv 

I. INTRODUCTION 

Cosmogenic neutrinos appear in any scenario proposed 
to explain the most energetic cosmic rays. In particular, 
if the observed air showers of up to 10 11 GeV Q are 
produced by primary protons, in their way to the Earth 
these protons will interact with the cosmic microwave 
background (CMB) photons and produce pions: 
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The flux of cosmogenic neutrinos would then be created 
in the decay of the charged pions, and it will appear 
correlated with observable fluxes of nucleons and photons 
(see @ for a recent review). 

Cosmogenic neutrinos are of great interest as probes of 
new TeV physics because they provide very large center 
of mass energies. In addition, the relative effect of new 
physics on the weakly interacting neutrinos is larger than 
on quarks or charged leptons, making it easier to see 
deviations. At these energies the new physics may be 
able to compete with the weak interactions and provide 
signatures that could be detected in deeply penetrating 
air showers and neutrino telescopes. 

In particular, one expects that at transplanckian en- 
ergies gravity dominates over all the other interactions. 
This will be the case when a cosmogenic neutrino inter- 
acts with a terrestrial nucleon in models with extra di- 
mensions and the fundamental scale Md at the TeV 0. 
The possibility of black hole (BH) formation p| by cos- 
mogenic neutrinos has been discussed by several groups 
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These analyses are based on a geo- 
metric cross section, which assumes gravitational collapse 
if the neutrino interacts at impact parameter distances 
smaller than the Schwarzschild radius Rs of the system. 
The collapse involves strongly coupled gravity and is not 
calculable perturbatively, but if Rs 3> M^ 1 (y/s Md) 
one expects that the estimate will not be off by any large 
factors ^(J. It is found, however, that the vN cross sec- 
tion is dominated by the low- a; region, with y/s at the 
parton level close to Md , and most of the BHs produced 
will have a radius Rs ~ Mp . In this regime the amount 
of gravitational radiation emitted during the collapse or 
the topology of the singularity are important effects that 
add uncertainty to the geometric estimate. 

Here we study the gravitational interaction at larger 
distances, where it can be calculated using the eikonal 
approximation [Til If3 . 1 1.1 IT^ . This approximation in- 
volves linearized gravity and is not affected by the un- 
certainties in the cross section for BH formation. In the 
next Section we show that, for the typical energy E v of 
cosmogenic neutrinos, the eikonalized vN cross section 
depends very mildly on how the theory is completed in 
the ultraviolet (UV), at energies around Md- We also 
show that the amount of gravitational radiation emitted 
during the scattering is small. In these processes the neu- 
trino interacts at impact parameters larger than Rs and 
therefore with a larger cross section than for BH pro- 
duction. At the parton level the neutrino scatters elas- 
tically and transfers a small fraction y of its energy to 
a quark or a gluon, which starts then a hadronic shower 
of energy yE v . Three main features characterize these 
processes and distinguish them from standard model or 
BH events. First, the shower has a typical energy much 
smaller than the energy (10 8 to 10 11 GeV) of the incom- 
ing neutrino. Second, a charged lepton is never produced 
in the starting point of the shower. Finally, the neutrino 
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is not destroyed in the interaction, it keeps going with 
essentially the same energy and may interact again. We 
show in Section III that neutrino telescopes are then ideal 
experiments to observe these elastic processes: they are 
designed to detect hadronic showers of energy down to 
100 TeV < 10 9 GeV (below 100 TeV the atmospheric 
background dominates), and their big volume (1 km 3 in 
IceCube |15|1 would favor multiple-bang events. In the 
final section we discuss and summarize our results. Our 
analysis here completes our work in 12], where aspects 
like the UV dependence of the cikonal amplitude, gravi- 
tational bremsstrahlung, or multiple-bang events in neu- 
trino telescopes were not discussed. 



II. TEV GRAVITY 

The simplest picture of TeV gravity includes only two 
free parameters: the value of the higher-dimensional 
Planck scale Md, and the number n of compact di- 
mensions where gravity propagates. A third param- 
eter, the (common) length 2irR of the n dimensions, 
could be deduced from the 4-dimensional Newton con- 
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At processes below Md the model-independent sig- 
nature of extra dimensions is graviton emission. The 
amount of energy radiated would be proportional to the 
accessible phase space or, in the Kaluza-Klein (KK) pic- 
ture, to the number of KK modes of mass below the 
center of mass energy. In this type of experiments for 
a given n one sets bounds on R and then deduces the 
limits on Md- From collider experiments one obtains 
M D > 1.4 (1.0) TeV for n = 2 (> 3) 16], whereas from 
SN1987A the bounds go up to 22 TeV for n = 2 [HJ. One 
should keep in mind, however, that the gravitons emit- 
ted in the supernova explosion have a KK mass below 
« 50 MeV. The simple picture with two extra flat di- 
mensions could be modified above 50 MeV, for example, 
with four more dimensions at R' ~ (100 GeV) -1 , which 
would bring the fundamental scale of gravity down to 1 
TeV without affecting the physics in the supernova. It 
could also be that some other mechanism (a warp fac- 
tor in |l8() gives an extra mass of order > 50 MeV to 
the KK excitations, invalidating all the bounds based on 
supernovas. 

The bounds obtained from transplanckian collisions 
are complementary in the sense that given n they are 
a direct probe of Md, and R is then adjusted in order 
to reproduce Gn- At energies above Md and impact 
parameters smaller than R the collision is a pure higher- 
dimensional process independent of the compactification 
details that fix the value of the effective Newton constant. 
The transplanckian collision does not see that the extra 
dimensions are compact, they could be taken infinite with 
no effect on the cross section. 



A. Neutrino-parton amplitude 

The TeV gravity model should be embedded in a string 
theory, which would relate Md with the string scale Mg. 
In the simplest set-up 19] the standard model (SM) 
fields (open strings) would be attached to a 4-dimensional 
brane, whereas gravity (closed strings) would propagate 
in the whole Z?-dimensional space. In this case 
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with g the string coupling. The transplanckian regime 
corresponds then to energies above the string scale, where 
any tree-level amplitude becomes very soft. In the ultra- 
violet string amplitudes go to zero exponentially at fixed 
angle and, basically, only the forward (long-distance) 
contribution of the graviton survives (the forward con- 
tribution of the SM gauge bosons also survives, but it 
is subleading above Md due to the smaller spin of the 
vector bosons). This is precisely the regime where the 
eikonal approximation is valid. 

Let us consider the elastic collision of a neutrino and a 
parton that exchange D-dimensional gravitons (see flU 
Il4| for details). The eikonal amplitude .4eik(s,i) resums 
the infinite set of ladder and cross-ladder diagrams. It is 
reliable as far as the momentum carried by the gravitons 
is smaller than the center of mass energy or, in terms 
of the fraction of energy y — (E v — E' v )/E v lost by the 
incoming neutrino, if y = —t/s «1 (s and t refer to the 
Mandelstam parameters at the parton level) . In this limit 
the amplitude is independent of the spin of the colliding 
particles. Essentially, _4 C ik is the exponentiation of the 
Born amplitude in impact parameter space: 



Aik(M) = y J d 2 be^ h (« 
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where x(s, is the eikonal phase and b spans the (bidi- 
mensional) impact parameter space. The Born am- 
plitude corresponds to ^4 e ;k(s,i) in the limit of small 
x(s, b) and, therefore, the eikonal phase can be deduced 
from the Fourier transform to impact parameter space of 

-4.Born(M) : 
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Our Born amplitude comes from the i-channel exchange 
of a higher-dimensional graviton: 
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where the integral over momentum qt along the extra 
dimensions (equivalent to the sum over KK modes) gives 
an UV divergence if n > 2. The magic of the eikonal am- 
plitude is that it will be well defined despite we obtain 
it from an UV dependent Born amplitude. To under- 
stand that, let us first evaluate x(s,&) using dimensional 
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FIG. 1: The real (dashed) and imaginary (dotted) parts and 
the moduli (solid) of the functions F n in the eikonal amplitude 
© for n = 2, 6. 



rcgularization. The Born amplitude becomes 
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which implies 
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Although the eikonal phase diverges at b = 0, the ampli- 
tude in Eq. (@J is insensitive to that: the contributions 
from the region b <C b c are quickly oscillating and tend 
to cancel. A e ik(s,q) can be written 



»4 e ik(s, q) = -iwsb 2 c F n (b c q) 



with 
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where q = y/— t, and the integration variable is v = b/b c 
(see Fig^). For q < b^ 1 this integral is dominated by 
impact parameters around b c , and for q > b~ x by a saddle 
point at b s = b c (n/qb c )^^ l+1 \ 



B. Cutoff dependence of the amplitude 

Let us now use an UV cutoff A to regularize the Born 
amplitude in Eq. (jSJ. This may be more physical than 
dimensional regularization since _D-dimensional gravity 
must be completed (embedded in a consistent theory) at 
energies around Md. For example, in the simple brane- 
world setting mentioned above all the KK excitations of 
the graviton with mass ( or q %) larger than the string scale 
decouple exponentially 20], an effect that would mimic 
a cutoff. In any case, we can use the cutoff to estimate 
the UV dependence of these eikonalized processes. 



The Born amplitude becomes 
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where I n (A/q) diverges like (A/q) n ~ 2 with the cutoff (the 
divergence is logarithmic for n = 2). It is now straight- 
forward to find 
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where -Ko(£) is a modified Bessel function of the second 
kind. Expanding l|12|l in powers of Ab we obtain 
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with A n {Ab) approaching the constant 2 2_,l /r 2 (n/2) for 
large values of Ab. This expression tells us that the cutoff 
introduces corrections to the eikonal phase which are rel- 
evant only at impact parameters b < A . This region in 
b gives a negligible contribution to ^4 e ik(s, q) (see Section 

liTr31 



C. Non-linear corrections and soft graviton 
emission 

The eikonal amplitude in Eq. 10 is well denned for all 
values of s and q. However, as q (or y = q 2 / s) grows non- 
linear corrections (H diagrams) become important |14| . 
The relevance of H diagrams implies a regime with strong 
gravitational coupling and important graviton emission 
(soft bremsstrahlung) . The strong coupling can be ex- 
pected just by inspecting the eikonal amplitude, since for 
—t/sm 1 the saddle point b s that dominates the integral 
in impact parameter space approaches the Schwarzschild 
radius Rs of the system: 
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A process with typical impact parameter b < Rs will not 
be properly described by the eikonal amplitude, since 
nonlinear corrections will be of order one. On the other 
hand, in eikonal processes with j « 1 the main contri- 
bution to the amplitude in Eq. Q comes from impact 
parameters much larger than Rs, where nonlinear effects 
are small (see Section III D(l . 

Soft graviton emission is also a consequence of nonlin- 
ear couplings, it appears as an imaginary contribution 
to the eikonal phase corrected by H diagrams (xh) plf . 
This contribution is of absorptive type, it damps the elas- 
tic cross section showing that a Bloch-Nordsieck mech- 
anism is at work. For a given value of b, the average 
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number N so f t of gravitons radiated during the scattering 
can be read directly from xh 13 H3 : 

(, \ 3n+2 
f J ' (15) 

where 

b r = (b?R 2 s n+2 )^ « (G^s 2 )^ 5 . (16) 

Therefore, the typical (transverse) momentum radiated 
will be Q pa A/s ft& _1 . Notice that to obtain the energy 
lost by the incoming neutrino this momentum must be 
boosted to the nucleon rest frame. In an eikonal scatter- 
ing the dominant impact parameter distance is (b) ~ b s . 
Both Q ps b- 1 ps M D ( 2 /M| ) /s) 1 / 2n + 2 and the number of 
gravitons N soit ps 2/ (3«+2)/(2„+2) (s / M 2 )(n +2)/(2„+2) de _ 

crease for decreasing values of y, implying that for y <C 1 
the amount of gravitational radiation during the scatter- 
ing is small (see below a numerical example). 

On the other hand, in a collision at (b) ss Rs one ex- 
pects a large fraction of energy transferred from the neu- 
trino to the parton, a large scattering angle, and a sig- 
nificant fraction of energy lost to radiation. At these and 
smaller values of b one would also expect black hole (BH) 
formation |j, la, ISl tl la l2l III! LUj< It has been shown, 
however, that a number of factors (angular momentum, 
charge, geometry of the trapped surface, radiation before 
the collapse) make a precise estimate difficult, specially 
for light BHs of mass just above Md- 

D. Numerical analysis of the uN eikonal cross 
section 

To understand the relative relevance of the different 
scales and processes involved, in this Section we will con- 
sider the scattering of a 10 10 GeV neutrino with a nu- 
cleon N = (n + p)/2. The vN center of mass energy 
is in this case y/s = s/lra^Ey = 141 TeV. We will take 
n = 2 or n = 6 extra dimensions and a fundamental scale 
Md = 1 TeV. The transplanckian regime will then in- 
clude the partonic processes of energy y/l = y/xs > Md, 
i.e., x > 5 x 10~ 5 . To evaluate the cross sections we 
will use the CTEQ5 parton distribution functions (PDFs) 
[22| . which are available both for fortran and Mathemat- 
ica codes. We will base our analysis on the kinematical 
variable y = [E v — El)/E„, which fixes q 2 = ys and the 
dominant impact parameter distance (b) in the eikonal 
process ({b) ps b s if q > b^ 1 or (b) bs b c if q < b^ 1 ). 
We will evaluate the PDFs at this dominant distance 

We find (see Fig. EJ) that the differential cross section 
^f* = / dx xs nbj\F n (b c q)\ 2 V ffan) (17) 

grows as y decreases [ll|. F° r example, for n — 2 (6) it 
is a factor of 265 (62) larger at y = 10~ 3 than at y = 0.1. 
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FIG. 2: Differential cross section da^ /dy of a 10 10 GeV 
neutrino for n = 2, 6 and an UV cutoff A — > oo (solid) and 
A = Md (dashed). 



The small y region corresponds to long distance processes 
where the neutrino interacts with a parton and transfers 
only a small fraction of its energy. This region is less 
important for n — 6 than for n = 2 extra dimensions, 
since then gravity dilutes faster and becomes weaker at 
long distances. On the other hand, values of y close to 
1 mean shorter distance interactions. Using Eq. 10 we 
can evaluate the contribution to the cross section from 
different regions in impact parameter space. We obtain, 
for example, that for a vN process with y = 0.5 a 52% of 
the eikonal cross section comes from impact parameters 
b < Rs if n — 2 (or a 71% if n = 6). In these pro- 
cesses with y ~ 1 the eikonal amplitude will be corrected 
by nonlinear contributions of the same order. Therefore, 
we will use the eikonal amplitude to evaluate elastic pro- 
cesses with y < y max ~ 0.2 only. Our results will depend 
very mildly on the actual value of y ma x, since the bulk 
of the cross section comes from the small y region. For 
example, for n = 2 the eikonal vN cross section for pro- 
cesses with 10~ 6 < y < y max is = 1.97 x 10~ 2 mb 
if y max = 0.4 or - 1.91 x 10~ 2 mb if y max = 0.1. 
For n — 6 the values of the cross would change with y max 
from 7.5 x 10~ 3 mb to 6.2 x 10~ 3 mb. 

The numerical relevance of the UV cutoff (see Section 
III Ut is expressed in Fig. |3 where we plot der^ /dy for 
A — > oo (solid lines) and for A = Mjj (dashed line). 
In the later case KK modes of the graviton heavier than 
Md are decoupled. We find that the vN differential cross 
section changes less than a 10% for 10~ 5 < y < 10~ 2 , so 
the cutoff dependence of the eikonal amplitude in these 
processes is not important. 

Let us now consider the geometric cross section (at the 
parton level) ctbh = T.R5, with Rs given in Eq. (|14fl . ctbh 
includes all the processes (elastic and inelastic) at impact 
parameter distances b below Rs, and it can be used to 
estimate the rate of black hole formation. As explained 
before, eikonal scatterings of small y will be dominated 
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by values of b larger than Rs- Therefore, the overlapping 
between these soft eikonal processes and the processes in 
the (inclusive) geometrical cross section will be negligible. 
It is then justified to consider two types of transplanck- 
ian (s > M|,) processes: elastic (long-distance) processes 
where the neutrino transfers to the partons a small frac- 
tion y < y max of its energy and keeps going, and shorter 
distance (b < Rs) hard processes where the neutrino loses 
in the collision most of its energy, possibly collapsing into 
a BH. To estimate the relative frequency of these two pro- 
cesses when a 10 10 GeV neutrino scatters off a nucleon, 
we can compare the eikonal cross section with y in- 
tegrated between 1(T 5 and 0.2 with <rg$. For n = 2 (6) 
we obtain = 1.94 x 1(T 2 mb (6.88 x 1(T 2 mb) and 



crg£ = 9.82 x 1CT 4 mb (4.07 x lO" 3 mb) 



i.e. 



the neu- 



trino will have 12.5 (1-64) interactions in which it trans- 
fers to the nucleon between 100 TeV and 2 x 10 9 GeV of 
energy, per each short distance (black hole) interaction. 
The total energy lost by the neutrino in these 12.5 (1-64) 
interactions is 



E r ih- — 



1 



dyyE„^fl 
dy 



= 1.00 x 10 9 GeV (5.04 x 10 s GeV). (18) 

Finally, let us comment on the amount of gravitational 
energy radiated in these eikonalized scatterings. In a par- 
ton process of energy ^/xs and inelasticity y the energy 
lost to radiation in the c.o.m. frame is 
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which in the nucleon at rest frame is 
E la d(x,y) = E* ad 



E v 



(19) 



(20) 



The average energy lost through soft bremsstrahlung per 
each short distance gravitational interaction is then 
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We find that, for n = 2 (6), during the 12.5 (1-64) eikonal 
processes the 10 10 GeV neutrino radiates soft gravitons 
with a total energy of 1.56 x 10 9 GeV (4.90 x 10 8 GeV). 



III. SIGNALS AT NEUTRINO TELESCOPES 

The flux of cosmogenic neutrinos depends on the pro- 
duction rate of primary nucleons of energy around and 
above the GZK cutoff Eqzk- It will always appear corre- 
lated with proton and photon fluxes that should be con- 
sistent, respectively, with the number of ultrahigh energy 
events at AGASA and HiRes |l| and with the diffuse 7- 
ray background measured by EGRET |2^| . 
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FIG. 3: Left: cosmogenic neutrino fluxes referred in the text 
as higher (solid), lower (dashed) and minimal (dotted). We 
plot the fluxes for one flavor $ = <f> ve + <f>p t and assume that 
all flavors have the same frequency. Right: corresponding 
number of downward (0 < 8 Z < 7r/2) cosmogenic neutrinos of 
each flavor. 



We will base our analysis on the two neutrino fluxes 
described in |2| (solid and dashed lines in Fig. [3J| . The 
first one saturates the observations by EGRET, whereas 
for the second one the correlated flux of 7-rays contribute 
only a 20% to the data, with the nucleon flux normalized 
in both cases to AGASA/HiRcs. The higher flux pre- 
dicts 820 downward neutrinos of each flavor with energy 
between 10 8 GeV and 10 11 GeV per year and km 2 , ver- 
sus 370 for the lower one. We will also comment on the 
minimal flux described in |24| (dots in Fig.|3J), where the 
cosmic ray events above the GZK cutoff are assumed not 
to be protons |25| . The proton events around and below 
Eqzk imply then just 35 downward neutrinos of each fla- 
vor with energy between 10 s GeV and 10 11 GeV per year 
and km 2 . 

To observe a hadronic event inside the telescope, a cos- 
mogenic neutrino must first survive as it crosses the at- 
mosphere and the ice (or water) above the detector. Its 
typical interaction length in a medium of density p is 
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(22) 



where a uN = <7§m + ct bh 1S * ne total cross section to have 
an interaction that destroys the neutrino. It is usual to 
express the length in terms of its depth: xq — pLo (i.e. 
one meter of water has a depth of 100 g/cm 2 ). Notice 
that we include in a vN both the SM and the short dis- 
tance gravitational interactions, but we ignore the soft 
(eikonalized) gravitational interactions because they take 
from the neutrino just a small fraction of its energy (the 
distortion in the flux of neutrinos that reach the detector 
that these interactions produce is negligible). 

A neutrino from a zenith angle 9 Z must cross a column 
density of material 



x{8 z 



dl p(l,9 z ) 



(23) 
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In practice, the path in the atmosphere is negligible and 
x(9 z ) is just the depth of the water or ice above the de- 
tector. The probability that it does not interact before 
reaching the detector is then 

P snTV (E„,9 z ) = e- x / x ° . (24) 

Once in the detector, the probability of an event is 

Pu*(E v )Kl-e- I " ,NA <! , (25) 

where L is the linear dimension of the detector and 
the total cross section. Therefore, the total number of 
events in the telescope in an observation time T is 

N = 2nAT fdEvYsJ^r f dcosd z P SUIV P int , (26) 

where A is the detector's cross sectional area and <j> Vi 
the neutrino flux. Before a complete numerical analysis 
we would like to discuss the possibility of multiple-bang 
events inside the detector. 



A. Multiple-bang events 

If L is similar or larger than the interaction length Lq, 
then the neutrino may interact more than once inside 
the detector. This is possible because the neutrino is not 
destroyed in the first eikonal interaction, it keeps going 
with basically the same energy and can interact again. 

Let us assume that w <r^ ^> crgjj, <tI<m an d let us 
neglect the amount of energy lost by the neutrino in each 
interaction. It is straightforward to find the probability 
of having exactly N interactions (bangs) in a length L: 

P N (L) = e- L ^ { -^^. (27) 

For example, the probability of having only one interac- 
tion would be P\{L) — exp (—L/Lq)(L/L ); for L<Lo 
we have P\{L) as L/Lq, but for L ^S> Lq this amplitude 
goes to zero (it is very unlikely to have only one interac- 
tion). Given L, the most probable number of interactions 
is N = L/ Lq, which is also the average number of inter- 
actions: 

OO j 

n=i L ° 

Notice that the probability of having any type of event, 
i.e., at least one interaction, is (see Eq. 

OO 

P(L) = J2 P N = 1 ~ e ~ L/L ° > ( 29 ) 

N=l 

whereas the probability that this event includes more 
than one interaction (a multiple-bang event) would be 
P{L)-P X {L): 

fmuit(i) = 1 - e~ L / L "(l + Lj Lq) . (30) 



Double-bang events could also be produced by SM or 
BH interactions. Within the SM, the second bang would 
correspond to the decay of the tau created in the first 
interaction. The probability that this happens would be 
the probability that a v T has a SM charged current inter- 
action times the probability that the tau lepton decays 
inside the detector. If the vN interaction results into a 
BH, its evaporation would also produce taus that could 
decay inside the detector. For the double-bang tau event 
to be contained inside a detector like IceCube (1 km of 
length with 125 m between strings), the energy of the tau 
lepton must be between 2.5 x 10 6 GeV and 10 7 GeV 0- 



B. Numerical example 

Let us consider again a single neutrino of E v = 10 10 
GeV. Within the SM, its interaction length in ice is 
Lq = 440 km. This means that typically it could reach 
the center of AMANDA or IceCube (1.8 km below the 
antarctic ice ^1|) from angles cos6* 2 > —0.03. If there 
are n = 2 (6) extra dimensions and Mjj = 1 TeV, the 
interaction length before a hard gravitational interaction 
would be just Lq = 17 km (4 km), which corresponds to 
cos6 z > 0.11 (0.44). For M D = 2.8 (4.5) TeV cr^ as 
and L as Lq M . 

If the neutrino reaches the detector, within the SM the 
probability that in L = 1 km (the longitudinal dimension 
of IceCube) it starts a hadronic shower is P-^f — 2.2 x 
10~ 3 . For a v r neutrino, the probability of an event with 
a tau lepton in the initial point of the hadronic shower is 
1.6 x 10~ 3 . 

If Mr> = 1 TeV the probability of a short distance 
gravitational event would be pPf* = 0.06 if n = 2 (or 0.22 
for n = 6). To find the probability of a soft eikonalized 
interaction we need to evaluate 

rVm*x drr L ' N 

•tf-L*^- (31) 

with y max = 0.2 and y min = (100 TeV)/E v (the energy 
of the shower should be above 100 TeV to avoid the 
atmospheric background). The probability of an event 
in a length L would then correspond to Oi nt w in 
Eq. (|25|) . which for L = 1 km and n = 2 (6) gives 
Jj" t k = 0.56 (0.33). This probability includes events 
with one bang: P" k = 0.36 (0.27), with two bangs: 
P| lk = 0.15 (0.06), and with more than two bangs: 
= 0.05 (0.008). 

As explained above, double-bang events could also be 
produced by SM interactions. The probability that a 
10 10 GeV v r produces a tau lepton of energy between 
2.5 x 10 6 and 10 7 GeV is just P 2 SM w 6.8 x 10" 5 . 

In Fig. 0] we express the probability P SUIV that the 10 10 
GeV neutrino survives through 1.8 km of ice to reach 
vertically IceCube for different values of Mr> (for large 
values of M D P SUTV « Pj^£, ~ 1). We also plot the prob- 
ability that if it has reached the detector it experiences 
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FIG. 4: The different probabiiities defined in the text for a 10 GeV neutrino reaching IceCube from 6 Z = as a function of 
Md for n — 2 and n — 6. 



a /iard interaction (P;™) or a so/f observable interac- 
tion (i^„ t k ) within a distance of L = 1 km. The prod- 
uct Povont = PsurvPnt would give the probability that 
the neutrino gives a signal in IceCube. We obtain that 
this probability is larger for eikonal than for BH events. 
We also find that it is maximal for Mo 0.8 TeV. For 
lower values of Mo the neutrino tends to interact before 
reaching the detector, and for larger values it tends to go 
through the detector without interactions. In Fig. 0] we 
also plot the probability of an eikonal event that includes 
only one bang (P^ vcnt ), two bangs (P e 2 vent ), or more than 
two bangs (P e >e„t)- If Md £ 1.5 TeV the probability of 
more than one bang would be very small. 

Within a distance L the average energy lost by the neu- 
trino in eikonal interactions and radiated through gravi- 
tational bremsstrahlung would be, respectively, 



(P eik ) = LpN A 



dy 



(P rad ) = LpN A 



dx 



Mils 



dyE, 



rad 



7 eik 



dx dy 



(32) 



(33) 



with P ra( j given in Eq. (|2*U|) . For L = 1 km and n = 2 (6) 
the 10 10 GeV neutrino will lose (P cik ) = 6 x 10 7 (1.2 x 
10 8 ) GeV to hadrons and (P rad ) = 9.2 x 10 7 (1.2 x 10 8 ) 
GeV to gravitational radiation. This means that, as it 
propagates in the detector, the energy loss in these soft 
processes is negligible. In a typical interaction length Lq 
of a hard interaction (where the neutrino will lose most 
or all of its energy) we find that {E ci] ,)/E y = 0.10 (0.05) 
and (E TaA )/E v = 0.16 (0.05). 

Finally, let us study what is the typical energy of the 
hadronic shower started by the 10 10 GeV neutrino. We 
will increase Mo to 2 TeV (to avoid double-bang events) 
and take n = 2 (6). If the neutrino reaches the detector 
at IceCube, the probability that it starts a shower of en- 
ergy between 100 TeV and 0.2E„ is P int « l- e - LpNArT '^ , 




0.01 



0.001 



10 5 10 6 10 7 10 s 10 9 
E sh [GeV] 

FIG. 5: Probability that a 10 GeV neutrino that reaches 
IceCube starts an eikonal shower in each interval of energies 
for Md = 2 TeV and n = 2, 6. For example, the neutrino has 
a probability of 0.024 of starting a shower of energy between 
100 and 1000 TeV if n = 2. 

0.2. We obtain P int = 0.070 (0.045). Now, the probabil- 
ity Ap nt that this event has an energy between E\ and 
E 2 is 



AP int = 



P 



hit 



J cik 



dy^- 
dy 



(34) 



where cr^t ~ a oik Wltn V 



100 TeV/E„ and y n 



In Fig. we divide the interval from 100 TeV to 2 x 10 9 
GeV in 5 bins and express the probability that the shower 
energy is in each of the bins. We observe that the typical 
value is between 10 5 and 10 s GeV for n = 2 and between 
10 6 and 10 9 GeV for n = 6. 



C. Bounds from air showers 

The fact that the typical energy P s h = yE v of the 
shower started in these processes is much smaller than 
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the energy E v w 10 10 GeV of the incoming neutrino has 
implications in air shower experiments. In particular, the 
absence of deeply penetrating showers in AGASA and 
Fly's Eye could exclude vN cross sections between 0.01 
and 1 mb. Notice, however, that these experiments are 
only sensitive to showers of very large energy, with E s h 
around or above 10 9 GeV. Since the eikonal cross sec- 
tions that we are considering reach a large size only for 
low values of y, the typical showers that they produce 
would be invisible in AGASA and Fly's Eye. This is in 
contrast to processes like BH formation Jl or other pro- 
cesses of strongly interacting neutrinos |25j, where most 
of the energy of the initial neutrino would be transferred 
to the shower. 

A precise analysis of the bounds on Md from eikonal 
processes in deep air shower experiments can be found in 
[l2| . The limits obtained there, between 1 TeV for n = 2 
and 1.5 TeV for n = 6, are essentially the same as the 
ones from BH production in p|. 




10 6 10 8 10 10 10 12 10 6 10 8 10 10 10 12 
E sh [GeV] E sh [GeV] 



D. Cosmogenic neutrinos at AMANDA and 
IceCube 

Let us now study the total number of events at 
AMANDA (0.03 km 2 and a length of 700 m) and Ice- 
Cube (1 km 3 ) for the neutrino fluxes in Fig. [31 

In the SM, for the higher flux (910 downward cosmo- 
genic neutrinos of each flavor per year and km 2 ), we 
would expect 1.32 contained events per year in IceCube. 
Of those, 0.38 would come from a neutral current and 
0.94 from a charged current (one third of the events of 
each lepton flavor). The distribution of energy of these 
events is given in Fig.E|(we show it despite the low statis- 
tics just for a comparison purpose). Around 0.008 of the 
0.31 tau events would decay and give double bangs in- 
side the detector. For the lower flux (around 410 cos- 
mogenic neutrinos of each flavor per year and km 2 ) , we 
would expect just 0.50 SM events per year inside the de- 
tector, 0.12 of them containing a tau lepton, and just 
0.003 double-bang events. The numbers for AMANDA 
can be easily obtained just multiplying by a volume fac- 
tor Vam/Vic ~ 0.02, namely 0.03 SM events per year 
for the higher flux and 0.01 for the lower flux. 

In a scenario with n = 2 (6) extra dimensions, for the 
higher flux we obtain a signal above the SM background 
(1.32 contained events per year in IceCube) if Md < 
5.6 (4.9) TeV, whereas for the lower flux we have a signal 
above the 0.50 SM events if M D < 4.8 (4.5) TeV. 

The event rate at IceCube and AMANDA for different 
values of Md and a minimum energy of the shower of 
100 TeV is plotted in Fig. [7| We give the number of 
short distance (BH) and of soft (eikonal) events. We also 
include the number of double-bang eikonal events, which 
is significant for low values of Md . 

The energy distributions of contained hadronic showers 
in IceCube for both fluxes, M D = 2 TeV and n = 2 (6) 
are also shown in Fig. H3 There is a clear difference be- 



FIG. 6: Energy distribution (events per bin) of the eikonal 
(solid), BH (dashed) and SM (dotted) events in IceCube per 
year for the higher and the lower cosmogenic fluxes, Md = 
2 TeV and n = 2, 6. 



tween the energy distribution of eikonal and BH or SM 
events: while these have a shape similar to the cosmo- 
genic flux, eikonal events are typically of much lower en- 
ergies, specially for n = 2. 



IV. SUMMARY AND DISCUSSION 

Cosmogenic neutrinos interact with terrestrial nucle- 
ons at center of mass energies \Z2ra^E v 100 TeV, 
so they can be used as probes of new TeV physics in 
neutrino telescopes. In particular, the possibility of BH 
formation in models with extra dimensions has been en- 
tertained by several groups. These analyses are based 
on a geometric cross section that assumes single BH pro- 
duction whenever the neutrino and the parton interact at 
impact parameters smaller than i?g. The problem with 
this estimate is that, despite the large energy of cosmo- 
genic neutrinos, the vN cross section is dominated by 
the small x region and most of the BHs produced in a 
neutrino telescope would be very light, with masses just 
above Mb- These light BHs would be very sensitive to 
effects like graviton emission during the collapse or non- 
thermal effects in the evaporation, which add uncertainty 
to any estimate. 

In this paper we have analyzed in detail a different type 
of signal. It is produced when the neutrino interacts elas- 
tically with a parton at typical distances larger than Rs 
and transfers a small fraction y of its energy. The pro- 
cess is properly described by the eikonal approximation. 
We have shown that the cutoff dependence of the eikonal 
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FIG. 7: Contained events per year in IceCube and AMANDA 
for the higher (thick) and the lower (thin) cosmogenic fluxes 
and n — 2, 6. We show eikonal (solid), multi-bang (dashed- 
dotted) and BH (dashed) events. 



cross section is small (see Fig. |2J), and that non-linear 
corrections and graviton emission are negligible. The dis- 
tinct experimental signature of these processes would be 
a hadronic shower of energy yE v <C E v . A muon, a tau, 
or an electromagnetic shower would never be produced 
in the initial vN interaction. 

If IceCube observes contained showers above 100 TeV, 
their energy distribution and the absence of charged lep- 
tons in the starting point of the shower would then suffice 
to decide whether they are due to SM or TeV-gravity in- 
teractions, as we argue below. 

Let us first suppose that the cosmogenic neutrino flux 
(important at energies between 10 8 and 10 11 GeV) is 
smaller than expected (see the minimal flux in Fig- EJ - 
In this case there will be no SM contained showers in 
IceCube. However, for sufficiently low values of Md (the 
optimal value is 0.8 TeV, see Pelent m Fig- El the eikonal 
cross section grows and even a very low cosmogenic flux 
could imply contained showers from gravitational inter- 
actions. These showers would be typically less energetic 
than the initial neutrino that produce them (see distri- 
bution in Fig. EJ. One may then wonder if these events 
could be due to a large flux of neutrinos in the interme- 
diate energy range (10 5 — 10 s GeV). The observation of 
an initial muon in a 24% of the showers would be char- 
acteristic of SM interactions in this range of energies, 
whereas the absence of muon events would be consistent 
with soft gravitational interactions of cosmogenic neutri- 
nos of much higher energy. In addition, the SM contained 
showers would always come together with a large number 
of muon events of similar energy produced outside the de- 



tector, which would be absent for TeV-gravity contained 
showers. 

Let us now suppose a flux of cosmogenic neutrinos 
within the expected limits (fluxes higher and lower in 
Fig- GIJ)- In this case IceCube will observe SM contained 
showers in the range 10 s to 10 11 GeV. If M D w 5 TeV the 
number of eikonal processes will be similar to the number 
of SM events. However, their energy distribution will be 
different and both types of processes can be separated. 
Again, the absence of charged leptons in the initial in- 
teraction point would distinguish these events from SM 
events of lower energies. For the higher flux in Fig. 
we expect more than one contained non-standard shower 
per year in IceCube if Md < 6.0 (5.5) GeV for n = 2 (6). 

We then conclude that for any intermediate-energy and 
cosmogenic neutrino fluxes, an enhanced rate of neutral 
versus charged current events of energies above 100 TeV 
could be explained by TeV-gravity interactions. 

These interactions could also produce a very peculiar 
signal for relatively low values of Md- In a typical eikonal 
process the neutrino loses a small fraction y of its energy 
and keeps going, so it can interact several times inside 
the detector. This effect is specially important for low 
values of n, where gravity dilutes slowly with the dis- 
tance (notice that for n < 2 it becomes a long-distance 
interaction and the total eikonal cross section is diver- 
gent). If n = 2 and Md < 0.9 TeV the average interac- 
tion length of a 10 10 GeV neutrino between two eikonal 
interactions of E > 100 TeV becomes smaller than the 
longitudinal dimension of IceCube. Therefore, we would 
expect two (or more) bangs of 10 5 to 10 8 GeV inside 
the detector. We think this type of events could be eas- 
ily distinguished from possible double-bang SM events, 
where the first bang corresponds to a v T N charged cur- 
rent interaction and the second one to a tau decay. First 
of all, the typical energy of the SM double-bang event 
would be necessarily between 2.5 x 10 6 and 10 7 GeV; for 
lower energies the tau decays before 125 m (the separa- 
tion between strings at IceCube) and for larger energies 
the second bang would be out of the detector. Second, 
there should be a clear trace in the detector as the tau 
propagates between the two bangs, which is absent in 
TeV-gravity events. 

In summary, we think that elastic eikonalized interac- 
tions provide a clear (distinguishable from possible SM 
events) and model-independent (insensitive to how the 
theory is completed in the UV) signal of TeV gravity. Be- 
ing at impact parameter distances larger than R$, these 
interactions have a cross section that is larger than the 
geometric cross section to produce a BH. The eikonal 
event would be much less energetic than a SM or a BH 
event, but neutrino telescopes are sensitive to showers of 
energy up to four orders of magnitude below the aver- 
age energy of cosmogenic neutrinos. The values of the 
fundamental scale of gravity that IceCube could reach, 
around 5 TeV, are comparable to those to be explored at 
the LHC or the ILC [l|. 

This work has been supported by MCYT (FPA2003- 



10 



09298-C02-01) and Junta de Andaluci'a (FQM-101). 
J.I.I, and D.M. acknowledge financial support from the 
European Community's Human Potential Programme 
HPRN-CT-2000-00149. We thank Eduardo Battaner, 



Tommaso Chiarusi, Francis Halzen, Marek Kowalski, 
Paolo Lipari, Teresa Montaruli, Sergio Navas, Andreas 
Ringwald and Christian Spiering for useful discussions. 



[1] L. Anchordoqui, T. Paul, S. Reucroft and J. Swain, Int. 

J. Mod. Phys. A 18 (2003) 2229. 
[2] D. V. Semikoz and G. Sigl, JCAP 0404 (2004) 003. 
[3] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. 

Lett. B 429 (1998) 263; I. Antoniadis, N. Arkani-Hamed, 

S. Dimopoulos and G. Dvali, Phys. Lett. B 436 (1998) 

257. 

[4] P. C. Argyres, S. Dimopoulos and J. March- 
Russell, Phys. Lett. B 441 (1998) 96; R. Emparan, 
G. T. Horowitz and R. C. Myers, Phys. Rev. Lett. 85 
(2000) 499; D. M. Eardley and S. B. Giddings, Phys. 
Rev. D 66 (2002) 044011. S. Dimopoulos and G. Lands- 
berg, Phys. Rev. Lett. 87 (2001) 161602; S. B. Giddings 
and S. Thomas, Phys. Rev. D 65 (2002) 056010. 

[5] J. L. Feng and A. D. Shapere, Phys. Rev. Lett. 88 (2002) 
021303; L. A. Anchordoqui, J. L. Feng, H. Goldberg and 
A. D. Shapere, Phys. Rev. D 65 (2002) 124027; Phys. 
Rev. D 66 (2002) 103002; Phys. Rev. D 68 (2003) 104025. 

[6] A. Ringwald and H. Tu, Phys. Lett. B 525 (2002) 135; 
M. Kowalski, A. Ringwald and H. Tu, Phys. Lett. B 529 
(2002) 1; S. I. Dutta, M. H. Reno and I. Sarcevic, Phys. 
Rev. D 66 (2002) 033002; A. Cafarella, C. Coriano and 
T. N. Tomaras, "Cosmic ray signals from mini black holes 
in models with extra dimensions: An analytical / Monte 
Carlo study" , |arXiv:hep-ph /0410358| 

[7] J. Alvarez-Muniz, J. L. Feng, F. Halzen, T. Han and 
D. Hooper, Phys. Rev. D 65 (2002) 124015. 

[8] E. J. Ahn, M. C avaglia and A. V. Olinto, 
|arXiv:hep-ph/0312249| 

[9] R. Casadio and B. Harms, Int. J. Mod. Phys. A 
17 (2002) 4635; D. Stojkovic, Phys. Rev. Lett. 94 
(2005) 011603; T. G. Rizzo, "Collider production 
of TeV scale black holes and higher-curvature grav- 
ity", |arXiv:hep- ph/0503163 J. L. Hewett, B. Lil- 
lie and T. G. Rizzo, "Black holes in many dimen- 
sions at the LHC: Testing critical string theory", 
arXiv:hep-ph/0503178 H. Yoshino and V. S. Rychkov, 
"Improved analysis of black hole formation in high- 
energy particle collisions", arXiv:hep-th/0503171 
[10] S. B. Giddings and V. S. Rychkov, Phys. Rev! D 70 
(2004) 104026. 

[11] R. Emparan, M. Masip and R. Rattazzi, Phys. Rev. D 
65 (2002) 064023; M. Masip, |arXiv:hep-ph/0210143 



[12] J. I. Illana, M. Masip and D. Meloni, Phys. Rev. Lett. 
93 (2004) 151102; D. Meloni, Acta Phys. Polon. B 35 
(2004) 2781. 

[13] G. »t Hooft, Phys. Lett. B 198 (1987) 61; I. J. Muzinich 
and M. Soldate, Phys. Rev. D 37 (1988) 359; D. Amati, 
M. Ciafaloni and G. Veneziano, Phys. Lett. B 197 (1987) 
81; D. Rabat and M. Ortiz, Nucl. Phys. B 388 (1992) 
570. 

[14] G. F. Giudice, R. Rattazzi and J. D. Wells, Nucl. Phys. 

B 630 (2002) 293. 
[15] J. Ahrens [IceCube Collaboration], 

arXiv:astro-ph/0305196 see also 

http://icecube.wis.edu/ 
[16] E A. Mirabelli, M Perelstein and M. E. Peskin, Phys. 

Rev. Lett. 82 (1999) 2236; G. F. Giudice, R. Rattazzi and 

J. D. Wells, Nucl. Phys. B 544 (1999) 3; for a review, see 

F. Feruglio, arXiv:hep-ph/0401033 
[17] S. Cullen andU. Perelstein, Phys. Rev. Lett. 83 (1999) 

268; S. Hannestad and G. G. Raffelt, Phys. Rev. D 67 

(2003) 125008 [Erratum- ibid. D 69 (2004) 029901]. 
[18] G. F. Giudice, T. Plehn and A. Strumia, Nucl. Phys. B 

706 (2005) 455. 
[19] S. Cullen, M. Perelstein and M. E. Peskin, Phys. Rev. D 

62 (2000) 055012; F. Cornet, J. I. Illana and M. Masip, 

Phys. Rev. Lett. 86 (2001) 4235. 
[20] I. Antoniadis and K. Benakli, Phys. Lett. B 326 (1994) 

69; S. A. Abel, M. Masip and J. Santiago, JHEP 0304 

(2003) 057. 

[21] D. Amati, M. Ciafaloni and G. Veneziano, Nucl. Phys. B 
347 (1990) 550. 

[22] H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 
12 (2000) 375. 

[23] P. Sreekumar et al, Astrophys. J. 494 (1998) 523. 

[24] Z. Fodor, S. D. Katz, A. Ringwald and H. Tu, JCAP 
0311 (2003) 015. 

[25] G. Domokos and S. Kovesi-Domokos, Phys. Rev. Lett. 82 
(1999) 1366; W. S. Burgett, G. Domokos and S. Kovesi- 
Domokos, "Low scale string unification and the highest 
energy cosmic rays", arXiv:hep-ph/0209162 Z. Fodor, 
S. D. Katz, A. Ringwald and H. Tu, "Strongly inter- 
acting neutrinos as the highest energy cosmic rays", 
|arXiv:hep-ph/0310iT2"l 



